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ABSTRACT 

Autonomous modal identification automates the 
calculation of natural vibration frequencies, 
damping, and mode shapes of a structure from 
experimental data. This technology complements 
damage detection techniques that use continuous or 
periodic monitoring of vibration characteristics. 
The approach shown in the paper incorporates the 
Eigensystem Realization Algorithm (ERA) as a 
data analysis engine and an autonomous supervisor 
to condense multiple estimates of modal parameters 
using ERA’S Consistent-Mode Indicator and 
correlation of mode shapes. The procedure was 
applied to free-decay responses of a Space Shuttle 
tail rudder and successfully identified the seven 
modes of the structure below 250 Hz. The final 
modal parameters are a condensed set of results for 
87 individual ERA cases requiring approximately 
five minutes of CPU time on a DEC Alpha 
computer. 

KEYWORDS: Modal testing, vibration testing, 
autonomous modal identification, damage 
identification using vibration testing, Eigensystem 
Realization Algorithm. 

INTRODUCTION 

Modal identification is the process of 
calculating natural frequencies, damping, and mode 
shapes (modal parameters) from experimental data. 
In most cases the experimental data are free-decay 
vibration measurements, frequency response 
functions (FRFs), or impulse response functions 
(IRFs) obtained by inverse Fourier transformation 
of FRFs or by using the Observer-Kalman Filter 
Identification (OKID) technique (Juang, 1994). The 
Eigensystem Realization Algorithm (ERA) 
identifies structural modal parameters from IRFs or 
free-decay data. It is a multiple-input/multiple- 
output, time-domain technique that efficiently 
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identifies many modes simultaneously including 
closely spaced and weakly excited modes. Since its 
introduction in 1984, ERA has been used 

effectively on a large number of structures, and 
many publications discuss both the theoretical and 
practical aspects of the method (Pappa, 1997a). 

ERA often generates good results using the 
default values of various algorithm parameters. 
However, depending on accuracy requirements, 
improvement of results may be necessary or 
desirable. One approach to improving modal 
identification results is to conduct additional tests. 
Weakly excited or weakly observed modes may be 
identified better by increasing excitation 
amplitudes, changing the number or orientation of 
exciters, improving signal-to-noise ratios of the 
data by different filtering, windowing, or averaging 
methods, adding or moving sensors, adjusting 
amplifier settings, etc. An entirely different 
approach is to perform additional ERA analyses of 
the same data set(s) using other values of the 
algorithm parameters. For example, the data record 
length may be increased to improve the accuracy of 
linear, lightly damped, low-frequency modes 
(probably at the expense of other modes with 
higher nonlinearity, damping, or frequency). The 
extent to which modal identification results 
fluctuate with changes in algorithm parameters 
depends on the dynamic complexity of the 
structure, notably its modal density and degree of 
nonlinearity, and on the quality of the experimental 
measurements. 

Autonomous modal identification of structures 
is a new technology in its early stages of 
development. In this paper, only the second 
approach for improving identification results is 
addressed; that is, the improvement of identified 
modal parameters by performing additional ERA 
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analyses of existing data set(s). In a previous 
study, autonomous modal identification was 
implemented using sliding-bandwidth filtering and 
both single-input and multiple-input analyses 
(Pappa et al, 1997b). In the present study, the size 
of the ERA Hankel data matrix and the assumed 
number of modes are varied. These are two 
additional, principal algorithm parameters. As more 
experience is gained with autonomous modal 
identification, a larger number of parameters will 
be included in the procedure. A full implementation 
will eventually incorporate all ERA algorithm 
parameters. 

In recent years, considerable attention has been 
given to the topic of damage detection of structures 
using vibration measurements (Doebling et al, 
1996; Flett, 1996; Grygier, 1994; James, 1996; 
Zimmerman et al, 1996). This technology applies 
to a wide range of structures including spacecraft, 
aircraft, buildings, bridges, off-shore oil platforms, 
and wind turbines. Many of the methods under 
development detect damage by monitoring changes 
in modal parameters. Using modal parameters for 
this purpose is attractive for several reasons, 
including 1) modal parameters change (to some 
degree) with any change in the mass, stiffness, or 
damping of the structure, 2) damage at inaccessible 
locations can be detected by changes in modal 
parameters measured at accessible locations, and 3) 
modal identification technology is widely used in 
laboratories and a broad experience base exists. The 
capability of generating structural modal 
parameters autonomously (i.e., without human 
involvement) for damage detection is the primary 
motivation of the work discussed in this paper. 

The next section of the paper describes the test 
article for this study, and the manner in which it 
was tested in the laboratory. The following section 
summarizes the autonomous modal identification 
technique. For this application, three sets of ERA 
analyses are performed comprising a total of 87 
cases followed by a final, overall mode 
condensation. The remainder of the paper shows 
principal identification results including a 
summary of the final, condensed set of modal 
parameters. 

TEST ARTICLE AND TEST DESCRIPTION 

Figure 1 shows the rudder/speed brakes (RSBs) 
of the Space Shuttle Orbiter. There are two pairs of 
RSBs, an upper pair and a lower pair. Each pair 
contains left and right surfaces that rotate in the 
same direction when used as a rudder and in 
opposite directions when used as a speed brake. The 
surfaces rotate up to plus or minus 27 degrees as a 
rudder, individually up to 49 degrees as a speed 
brake, or combined for joint rudder/speed brake 
control. The vertical tail structure is covered with 
reusable thermal insulation and is designed to 
withstand a 163 decibel acoustic environment with 


a maximum temperature of 350 degrees F. 

The structure used in this study is shown in 
Fig. 2. It is an upper-right RSB, and was tested 
with free-free boundary conditions by suspending it 
from the actuator attach points using springs and 
cables. The test occurred at the NASA Johnson 
Space Center (JSC) in December 1996. This 
structure is one component of a circa-1980 Vertical 
Stabilizer Assembly (VSA) being used at JSC and 
the University of Houston for finite-element model 
correlation and damage identification studies. The 
VSA originally served as an acoustic and static 
fatigue test article. The components of the VSA are 
being tested individually and assembled to develop 
techniques for correlating component and assembled 
models. .The structure will then be damaged in 
various ways to test damage identification methods. 
The RSB is approximately 48” x 7” x 120” in size, 
and weighs 1 80 lbs. 



Fig. 1. SPACE SHUTTLE ORBITER 




Fig. 2. MODAL TEST CONFIGURATION 
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Figure 3 shows the excitation and response y direction (perpendicular to the surface of the 

degrees-of-freedom used in the modal test. The structure) which emphasizes the primary, out-of- 

structure was excited several times at each of the 7 plane bending and torsional modes of the structure, 

locations shown in Fig. 3(a) with a force- In addition to calculating FRFs, the free-decay 

instrumented hammer, and frequency response vibration responses for the final impact at each 

functions (FRFs) were calculated between the excitation location were also recorded. The 

applied forces and each of the 47 accelerometers autonomous modal identification will analyze these 

shown in Fig. 3(b). Note that the majority of free-decay responses, 

impacts and response measurements are in the 



(a) 7 Impact Locations (b) 47 Accelerometers 

Fig. 3. EXCITATION AND RESPONSE LOCATIONS 



Time, sec Frequency, Hz 

(a) Excitation: 603y, Response: 672y 



(b) Excitation: 685y, Response: 61 1y 
Fig. 4. TYPICAL FREE-DECAY RESPONSES AND SPECTRA 
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Typical free-decay measurements and their 
corresponding FFT-magnitude spectra appear in 
Fig. 4. Each time history contains 2033 data 
points at a sampling rate of 800 samples/sec, and 
the data is alias-free out to a frequency of 312.5 
Hz. The relatively noisy data in Fig. 4(a) is 
representative of most measurements obtained 
using y-direction excitation at the actuator attach 
points (points 601-604). The relatively cleaner data 
in Fig. 4(b) is representative of most 
measurements obtained using y-direction excitation 
at the apex of the structure (point 685). Although 
not shown, most measurements obtained using x- 
and z-direction excitation at points 601 and 603, 
respectively, are even noisier than that shown in 
Fig. 4(a). The autonomous modal identification 
will analyze all 329 response measurements 
simultaneously to obtain global estimates of the 
modal parameters (329 responses = 47 

accelerometers x 7 impact locations). By using all 
data simultaneously in the analysis, the cleaner 
measurements tend to dominate in the calculations 
over the noisier measurements because of their 
larger amplitudes. 



Fig. 5. AVERAGE POWER SPECTRUM 

Figure 5 shows the average power spectrum of 
all 329 free-decay responses. It provides a good 


overview of the natural frequencies of the structure 
and their relative strengths in the measurements. 
This function is the average squared magnitude of 
all FFT spectra such as shown on the righthand 
sides of Figs. 4(a) and 4(b). It indicates that there 
are seven predominant modes of the structure from 
zero to 250 Hz (the frequency range of interest). 
There is also a small inflection in Fig. 5 at 72 Hz 
indicating an additional mode in this frequency 
range, but it is too poorly excited and/or too 
poorly sensed to be accurately identified. (This 
peak at 72 Hz is more apparent in other individual 
spectra not shown.) This mode is probably a local 
mode of the cable suspension system. 

AUTONOMOUS MODAL IDENTIFICATION 

Figure 6 is a flow diagram of the autonomous 
modal identification approach. There are four 
principal steps: 1) running the Eigensystem 

Realization Algorithm (ERA) (Pappa, 1997a), 2) 
applying thresholds to the results which filter out 
unreasonable values, 3) condensing (i.e., 
reconciling) the most-recent set of candidate modal 
parameters (the output of the thresholding step) 
with results obtained in previous passes around the 
feedback loop shown in Fig. 6, and 4) selecting 
appropriate algorithm parameters for the next ERA 
analysis. Note that the Mode Condensation block 
is also known as the Supervisor because it is 
responsible for determining how multiple ERA 
analyses are combined into a coherent, overall set 
of results. The following paragraphs discuss each 
of the four steps in more detail. The dashed line in 
Fig. 6 from the Supervisor to the Case Selector 
indicates a desirable feature that is under 
development, but it is not yet available. The final 
set of modal parameters is the last-computed 
output of the mode condensation process. 

ERA 

The Eigensystem Realization Algorithm serves 
as a data analysis “engine.” ERA is a well- 
established, documented procedure for estimating 
structural modal parameters (natural frequencies, 
modal damping, and mode shapes) from free-decay 
vibration responses (or impulse response 
functions) based on a set of user-specified 
algorithm parameters (Pappa, 1994). Although 



Fig. 6. AUTONOMOUS MODAL IDENTIFICATION APPROACH 
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default values and/or guidelines are available for 
each algorithm parameter, users have considerable 
freedom in selecting their values. With 
experimental data, the ERA results will fluctuate 
to some degree with any change in the algorithm 
parameters. The magnitude of fluctuation varies 
from mode to mode, and it depends in a complex 
manner on the characteristics and quality of the 
data, including its modal density, degree of 
nonlinearity, noise level, modal response 
amplitude, damping level, spatial linear 
independence of closely spaced modes, data 
sampling frequency, etc. 

In addition to generating modal parameters, 
ERA also calculates various “accuracy indicators” 
associated with each result. The principal ERA 
accuracy indicator is the Consistent-Mode Indicator 
(CMI) which provides a single, overall value for 
each mode ranging from zero to 100% (Pappa et al, 
1993). Modes with CMI values greater than 
approximately 80% are identified with high 

confidence. Modes with values from 80% to 1% 
display moderate to large uncertainty. Extraneous 
“noise modes” have CMI values of approximately 
zero. The CMI calculation uses both the temporal 
and spatial consistency of the results, and it is a 
robust, proven indicator of modal identification 
accuracy. 

Thresholding 

The thresholding step in Fig. 6 eliminates 
unfeasible modal identification results generated by 
ERA. A normal characteristic of ERA is that it 
generates extra, nonphysical modes in addition to 
the true, structural modes. The “noise modes” 


mentioned in the previous paragraph usually 
comprise a large fraction of these nonphysical 
modes. For good results with experimental data, at 
least half of the total number of identified modes 
will normally be nonphysical modes. Currently, 
nonphysical modes are eliminated from the results 
using the simple criteria shown in Fig. 7. Note 
that the CMI cutoff of 50% will allow some real, 
structural modes with poor accuracy to be deleted, 
but this is acceptable in most cases. (Hopefully, 
that mode is identified better in another ERA 
analysis.) The damping and frequency cutoffs of 
10% and 1%, respectively, are good selections for 
typical engineering structures having modal 
damping ratios on the order of 2% or less. These 
values would have to be increased for structures 
with higher damping. 


Thresholding 

Delete modes with: 

CMI < 50%, Damping < 0 or > 10%, 
or Frequency within 1 % of edges of 
analysis bandwidth. 


Fig. 7. THRESHOLDING LOGIC 

Mode Condensation (Supervisor) 

Figure 8 shows the mode condensation logic, 
which is presently the principal autonomous aspect 
of the procedure. Mode condensation selects the 
best estimate of each structural mode from 


Transfer Mode n from 



Fig. 8. MODE CONDENSATION LOGIC 


5 




multiple estimates of the mode. In other words, 
given a set of best estimates of the modal 
parameters (a set which is available initially after 
the first pass through the loop of Fig. 6), mode 
condensation decides whether or not new estimates 
of the modal parameters (referred to as Set 1 in 
Fig. 8) will be added to or replace results already 
in the best set (referred to as Set 2). 

Mode condensation consists of the following 
steps: 

1. For each mode in Set 1 (Mode n), compare 
its frequency with the frequencies of all modes in 
Set 2. If there are no modes of Set 2 within 10% 
in frequency of Mode n, add Mode n to Set 2. 
Recall that Mode n has already passed through the 
Thresholding step, so it is probably a true, 
structural mode (albeit with some inaccuracy if its 
CMI value is not high). 

2. For each mode of Set 2 within 10% in 
frequency of Mode n, calculate the Modal 
Assurance Criterion (MAC) (Pappa, 1994, p. 1 13). 
MAC is the square of the correlation coefficient of 
the two modal vectors (mode shapes) and measures 
the similarity of mode shapes. If all MAC values 
are less than 70%, add Mode n to Set 2 because its 
mode shape is significantly different than all 
modes in Set 2 with similar frequency. 

3. If there are one or more MAC values of 
70% or greater, select the mode in Set 2 having 
the highest MAC (Mode m), and compare its CMI 
value with the CMI of Mode n. If the CMI of 
Mode n (the new mode) is greater than that of 
Mode m, replace Mode m in Set 2 with Mode n. 

As mentioned earlier, information generated 
during the mode condensation activity will 
eventually be used to recommend algorithm 
parameters for subsequent ERA analyses (indicated 
by the dashed line in Fig. 6). A large amount of 
information is obtained during mode condensation 
concerning the characteristics of the data (such as 
the linear independence of closely spaced modes), 
and this information could be used to select better 
algorithm parameters. This capability is under 
development and is not yet available. 

Case Selector 

The final block in Fig. 6 is the Case Selector. 
Because the rudder/speed brake test article is a 
relatively simple structure, the case selection 
logic shown in Fig. 9 proved to be adequate for 
this application. Set A is a standard, default type 
of initial ERA analysis. Although it involves 35 
cases, the total execution time is quite short 
(approximately 75 secs) because the Hankel matrix 
and the number of measurements are fairly small. 
The final results of Set A (discussed in the next 
section of the paper) contained 2 modes with CMI 


values less than 90% (out of a total of 7 modes 
below 250 Hz). Accordingly, Set B was then 
executed in an attempt to improve the accuracy of 
these modes. After both Sets A and B completed. 
Set C combined the two groups of results into a 
final, overall set of results using the mode 
condensation logic of Fig. 8. 


Case Selection 

Set A (35 cases): 47 x 7 MEMO analysis 
with Hankel matrix of 470 rows x 70 cols. 
Assumed No. of Modes - 1 through 35. 

Stop here if all modes in the final results of 
Set A have CMI > = 90%. If not, proceed 
to Sets B and C. 

Set B (52 cases): 47 x 7 MEMO analysis 
with Hankel matrix of 470 rows x 105 cols. 
Assumed No. of Modes - 1 through 52. 

Set C: Overall Results - Combine results 
of Sets A and B by mode condensation. 


Fig. 9. CASE SELECTOR LOGIC 

Sets A and B are multiple-input/multiple- 
output (MIMO) analyses using all 329 free-decay 
response measurements simultaneously. Both sets 
use a Hankel data matrix with 470 rows (10 shifts 
of the 47 outputs). However, the number of 
columns in the matrix increases from 70 in Set A 
(10 shifts of the 7 inputs) up to 105 in Set B (15 
shifts of the 7 inputs). The maximum number of 
identified modes in any ERA analysis is always 
equal to one-half the minimum dimension of the 
Hankel matrix. For good results with experimental 
data, the maximum number of identified modes is 
typically in the range of approximately 2 to 10 
times the number of structural modes in the data 
bandwidth (up to 400 Hz in this problem). The 
numbers of columns in Sets A and B (70 and 105) 
were chosen based upon an estimated number of 
structural modes from zero to 400 Hz of 
approximately 15 obtained by counting the peaks 
in the average power spectrum. Fig. 5. 

MODAL IDENTIFICATION RESULTS 

Natural Frequencies and CMI 

Figures 10(a) and 10(b) show the ERA- 
identified natural frequencies between zero and 250 
Hz (the frequency range of interest) as a function 
of the assumed number of modes for Sets A and B, 
respectively. There are seven predominant modes 
of the structure in this frequency range as 
indicated. Each row of results in Fig. 10 

corresponds to a separate ERA analysis with the 
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(b) Set B (52 cases): Hankel Matrix Size - 470 rows x 105 columns 
Fig. 10. IDENTIFIED FREQUENCIES VS. ASSUMED NUMBER OF MODES 
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specified number of assumed modes. Each identified 
mode is represented by a short vertical dash at the 
associated frequency. The height and width of each 
dash express the confidence of the result based on 
the Consistent-Mode Indicator (CMI) (Pappa et al, 
1993), as described in the following paragraph. 
(Note: These natural frequency results occur prior 
to the Thresholding block of Fig. 6. For clarity, 
modes with CMI values less than 1% are not 
shown. ) 

The height of each vertical dash in Fig. 10 is 
directly proportional to the CMI value such that a 
mode with a value of 100% has a vertical dash 
height equal to the distance between y-axis tic 
marks. The width of each dash also increases with 
increasing CMI. This is done by using four line 
widths ranging from narrowest to widest for the 
following ranges of CMI: 1) below 80%, 2) from 
80% to 90%, 3) from 90% to 95%, and 4) above 
95%. (On a computer terminal, the four line widths 
correspond to four different colors, rather than line 
widths, to allow quick overview of the results.) 

In this application, the CMI accuracy indicator 
did an excellent job distinguishing the true, 
structural modes from the nonphysical, 
computational modes generated by ERA. Figures 
11(a) and 11(b) show typical results from Sets A 
and B, respectively. The seven structural modes 
have CMI values that are consistently much higher 
than those of the nonphysical modes. The majority 
of nonphysical modes, in fact, have values less 
than 1%. The unusually high values near 0 Hz, and 
at 120 and 240 Hz, are due to small DC offsets in 
the data and electrical noise, respectively. The 
other unusually high value at approximately 170 
Hz is probably a local mode of the cable 
suspension system. 

Although CMI clearly distinguishes the seven 
structural modes from the nonphysical modes, a 
close examination of Fig. 10 shows considerable 
variations of CMI for most of the structural modes 
as a function of the assumed number of modes 
(indicated by varying line widths). Out of the 14 
structural-mode results shown in Fig. 10(a,b), only 
Modes 2 and 3 in Fig. 10(b) have uniformly 
increasing CMI with increasing numbers of 
assumed modes. This desirable, convergent 
behavior occurs for wel 1 -exci ted, linear modes 
when the Hankel matrix size is “large enough.” On 
the other hand, all of the other results in Fig. 10 
display fluctuating CMI values as the assumed 
number of modes increases. Unfortunately, such 
fluctuations are the rule rather than the exception 
with many experimental data sets due to 
unavoidable “real world” effects, particularly 
nonlinearity and noise. The thresholding and mode 
condensation logic discussed earlier in the paper 
are designed to sift through the results of these 87 
cases and extract the most-confident set of 


structural modal parameters. 



F r eq uen c y , Hz 

(a) Set A: Top Row of Fig. 10(a) 



Fig. 11. TYPICAL DISTRIBUTIONS OF 
CMI VALUES 


Mode Condensation Results 

As indicated in Fig. 6, the final set of modal 
parameters from the autonomous modal 
identification approach is the output of the Mode 
Condensation process at the end of each case. If a 
sufficient collection of cases is run, the final, 
condensed set of modal parameters will be the best 
estimates of the true, structural modes (and 
nothing else). 

This autonomous modal identification approach 
performed extremely well in this application. 
Figure 12 ^shows the number of identified modes 
(i.e., the number of modes in the final, condensed 
set of results) at the end of each case for both Sets 
A and B. Recall from Fig. 9 that condensed results 
are obtained individually for Sets A and B, and 
then the individual results are combined in Set C 
using a final mode condensation. The dips in Fig. 
12 show that modes are being replaced according to 
the logic of Fig. 8. At each abscissa value, the 
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larger of the two numbers is the total number of 
modes in Set 2 at the completion of the case. The 
smaller of the two numbers is the number of modes 
retained from the previous case. In other words, the 
difference between the two numbers is the number 
of modes replaced as a result of the ERA analysis 
for the particular Case Number. The total number 
of identified modes in both Sets A and B steadily 
increases up to the true value of seven, with many 
fluctuations along the way indicating continuous 
improvement of results. This convergent behavior 
with continuous improvement of the results is 
precisely the way in which autonomous modal 
identification is designed to operate. 




Fig. 12. NUMBER OF IDENTIFIED MODES 
VS. CASE NUMBER 


Final Modal Identification Results 

The final set of results for both Sets A and B 
consisted of estimates of the seven true, structural 
modes. Note in Fig. 12 that both Sets converged to 
a stable number of modes (i.e., the replacement of 
modes slowed down and then ceased entirely as the 
Case Number approached 35 in Set A and 52 in Set 
B). Following completion of Sets A and B, Set C 
executed a final mode condensation on these two 
individual sets of results. 

Table 1 shows these final, overall modal 
identification results obtained in Set C. The second 
column in the table gives the Case No. from which 
the result was obtained. Note that 6 of the 7 modes 
came from Set B which used a larger Hankel data 
matrix size than Set A. In many cases, though not 
always, increasing the size of the ERA Hankel 
matrix will increase modal identification accuracy. 
However, there is a disproportionate increase in the 
computation time. For example, in this application 
Set A required 75 secs of CPU time on a DEC 
Alpha 4000 computer, whereas Set B required 213 
secs of CPU time but achieved higher accuracy in 
many of the modes. 

All of the CMI values in Table 1 are greater 
than 80% indicating high confidence in the results. 
Furthermore, five of the CMI values are above 95% 
indicating extremely high confidence in these 
modes. The righthand column in Table 1 also gives 
the Weighted Modal Phase Collinearity (MPC-W) 
which is a parameter used in the CMI calculation 
(Pappa, 1993). MPC-W measures the nearness of 
the mode shape to a monophase vector (i.e., to a 
classical, normal mode). Values greater than 99% 
are extremely high. 

Figure 13 shows the corresponding mode 
shapes. All seven of the modes are global, out-of- 
plane bending and torsional modes as expected. 
Each mode has a well-defined, predictable shape, 
substantiating the accuracy of the results. 
Furthermore, all seven frequencies and mode shapes 
agree closely with another, independent set 


Table 1. FINAL MODAL PARAMETERS AND ACCURACY INDICATORS (SET C) 


Mode 

Case 

Frequency, 

Damping 

CMI, 

MPC-W, 

No. 

No. 

Hz 

Factor, % 

% 

% 

1 

31 - Set B 

54.3 

0.17 

97.24 

99.94 

2 

44 - Set B 

69.3 

0.59 

98.52 

99.99 

3 

26 - Set B 

115.3 

0.35 

98.77 

99.94 

4 

45 - Set B 

158.8 

0.28 

96.13 

99.92 

5 

48 - Set B 

174.6 

0.73 

95.53 

99.58 

6 

42 - Set B 

226.6 

1.12 

87.61 

98.68 

7 

28 - Set A 

235.3 

0.61 

91.96 

95.89 
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(a) First Bending at 54.3 Hz 



(c) Second Bending at 115.3 Hz 



(e) Third Bending at 174.6 Hz 



(b) First Torsion at 69.3 Hz 



(d) Second Torsion at 158.8 Hz 



(f) Fourth Bending at 226.6 Hz 



(g) Third Torsion at 235.3 Hz 


Fig. 13. MODE SHAPES (SET C: FINAL RESULTS) 
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of results obtained by analyzing frequency response 
functions with the Polyreference technique. Good 
correlation of the two independent sets of results 
provides additional confidence in the accuracy of 
the autonomous modal identification approach 
described in this paper. 

CONCLUSIONS 

Considerable research is currently underway on 
several fronts in the area of damage detection of 
structures using changes in vibration 
characteristics. Many of the techniques under 
development require periodic calculation of 
structural modal parameters from experimental 
data. An attractive approach to this problem is 
autonomous modal identification. This new 

technology also promises improved speed and 
accuracy for more-traditional laboratory 

applications of modal identification techniques. 

The autonomous modal identification approach 
described in the paper is a fairly straightforward 
extension of the Eigensystem Realization 

Algorithm. It successfully identified the seven 
vibration modes of the Shuttle Orbiter tail rudder 
from zero to 250 Hz in approximately five minutes 
of computer time. As more experience is gained 
using this technology on a variety of other 
structures, additional logic and features will be 
added to the software to improve its performance 
and reliability in generating optimum modal 
parameters without human involvement. 
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